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ABSTRACT 
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Aims. Spiral structure is a prominent feature in many disk galaxies and is often outlined by bright, young objects. We study the 
?H distribution of young stellar clusters in grand-design spiral galaxies and thereby determine whether strong spiral perturbations can 

influence star formation. 

Methods. Deep, near-infrared JHK-maps were observed for ten nearby, grand-design, spiral galaxies using HAWK-I at the Very 
Large Telescope. Complete, magnitude-limited candidate lists of star-forming complexes were obtained by searching within the K- 

^ band maps. The properties of the complexes were derived from (H-K)-(J-H) diagrams including the identification of the youngest 

complexes (i.e. H Myr) and the estimation of their extinction. 

Results. Young stellar clusters with ages ^7 Myr have significant internal extinction in the range of Av=3-7'", while older ones 

O' ' typically have Av< 1™. The cluster limrinosity function (CLI^ is well-fitted by a power law with an exponent of around -2 and 

displays no evidence of a high luminosity cut-off. The brightest cluster complexes in the disk reach luminosities of Mr = -15?5 or 
estimated masses of 1O*M0. At radii with a strong, two-armed spiral pattern, the star formation rate in the arms is higher by a factor of 

• 2-5 than in the inter-arm regions. The CLF in the arms is also shifted towards brighter Mk by at least 0"4. We also detect clusters with 
i^H colors compatible with Large Magellanic Cloud intermediate age clusters and Milky Way globular clusters. The (J-K)-Mk diagram 

of several galaxies shows, for the brightest clusters, a clear separation between young clusters that are highly attenuated by dust and 
JL older ones with low extinction. 

^ Conclusions. The gap in the (J-K)-Mk diagrams implies that there has been a rapid expulsion of dust at an age around 7 Myr, possibly 

triggered by supemovae. Strong spiral perturbations concentrate the formation of clusters in the arm regions and shifts their CLF 
^ towards brighter magnitudes. 

Key words, galaxies: spiral - galaxies: structvu-e - galaxies: star clusters: general - galaxies: star formation - infrared: galaxies - 
I techniques: photometry 
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<3\ 1. Introduction tivestrengthof shocks in arms. This led them to suggest that den- 

sity waves trigger star formation in the vicinity of spiral arms. 

*0 A large fraction of the star formation in the local universe occurs tt, » ^ ..u a ■ ^ * 

• . . , ,. , , . . ,. ,, , , • , The star formation m the arm and mter-arm regions of three 
in gas-rich, disk galaxies as indicated by the Schmidt-Kennicutt • , j- j u i * i /omm • 3 * ^ ™ tt, 

^ , • ,%c,r^ ni J i^^r^x .1 , ■ r spirals was studicd Dy Foylc et al. (2010) usiHg data from thc 

law (Kennicutt 1989; Schmidt 1959), which relates the star tor- i» • i , , j- Tn, i j »u » . ^ ^ ■ ■ ^ 

^r; .or,„N , , , , • , ultraviolet to radio. They found that star formation in the inter- 

Cn mation rate (SFR) to the surtace/volume density ot gas at kpc vc .ji n u j-.u 

^ ^ , ™ ^ „ J- • , , , , , r^T-T. ^ arm regions was significant and only marginally enhanced in the 

^— I scales. The mnuence ol spiral structure on the global SFR ol ^^^^ 

galaxies is limited as indicated by the small difference seen 
.!h between multi-armed and grand-design spirals (Elmegreen & A study of cluster populations in 21 nearby spiral galaxies 
bx^ Elmegreen 1986). On smaller scales, the relation is less clear, led Larsen & Richtler (1999) to conclude that the age distribu- 
H partly owing to the low spatial resolution of gas maps (Bigiel tio" of clusters in disk galaxies shows no obvious peak, indi- 
et al. 2008), and may be influenced by local features such as res- eating fliat massive clusters are formed as part of an ongoing 
onances, shear, and spiral structures. This may indicate that the process rather than in bursts. They also concluded fliat the ra- 
main eff^ect of spiral structure is to concentrate the SFR into the dial distribution of clusters follows that of the H„ surface bright- 
arms rather than to change the global rate (Kennicutt 1998). "^ss. Hubble Space Telescope (HST) high resolution observa- 

„ . „ , , J , ni-T, ^ -.r, • i i tious (Larseu 2009) suggcstcd that thc initial clustcr mass fuuC" 

Seigar & James (2002) analyzed the SFR ot 20 spiral galax- ^. . ^ , . , ,. , , j i j o i, u. 

^ , J. , , r,r-n ■ 1 J • c 1 • .1 tiou m preseut-day spiral disks can be modeled as a Schechter 

les, concluding that the SFR is enhanced significantly in the ^ u u/ im/c\ -tu * fr™ o mSn/r 

... ^ , ? , J J , ^ „ • , , , function (Schechter 1976) with a cut-off mass Me » 2x10 Mp). 

vicmity ot their K-band arms and correlates well with the rela- ^ 

There are only a small number of grand-design spirals with 

prominent arms at a distance where individual clusters can 

Send ofprint requests to: P. GTOsh0\ be analyzed. One such galaxy, NGC 5194, was studied by 

* Based on observations collected at the European Southern Scheepmaker et al. (2007, 2009). These authors found that the 
Observatory, Chile; program: ESO 82.B-0331 spatial distribution of the star clusters younger than 10 Myr 
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Table 1. List of galaxies. Name and Hubble type are listed as given in RSA. Filter bands observed, seeing in the final, stacked K-band image 
expressed as the full-width at half-maximum (fwhm) in arcsec, and Ks surface brightness fis in magnitudes per pixel for which a signal-to-noise 
ratio of five was reached are also shown. The adopted distance modulus, m-M, distance D in Mpc, and the corresponding linear resolution As in 
pc are provided. The position angle PA and inclination angle LA adopted for the galaxies are listed in degrees. Limiting magnitudes K;, H;, J;, and 
Y; are indicated corresponding to a completeness of around 90% (see text). Finally, the total number of sources Nj in the galaxy disks for which 
aperture magnitudes could be estimated is given. 



Galaxy 


Type 


Bands 


fwhm 




m-M 


D 


As 


PA 


lA 


K, 


H, 


J/ 


Y, 


N, 


NGC 157 


Sc(s) I-II 


YJHK 


0.37 


23.9 


31.28 


18.0 


32 


40" 


45" 


20.19 


20.71 


21.27 


21.75 


2254 


NGC 1232 


Sc(rs) I 


YJHK 


0.38 


24.1 


31.48 


19.8 


36 


90* 


30* 


20.62 


20.80 


21.72 


22.21 


3177 


NGC 1300 


SBb(s) L2 


YJHK 


0.55 


24.0 


31.46 


19.6 


52 


87^ 


35"^ 


19.80 


20.30 


21.18 


21.30 


614 


NGC 1365 


SBb(s) I 


YJHK 


0.40 


23.7 


31.62 


21.1 


41 


40'' 


40'' 


20.17 


20.32 


21.24 


21.73 


2417 


NGC 1566 


Sc(s) I 


YJHK 


0.44 


23.1 


31.55 


20.5 


44 


41" 


27'' 


19.29 


19.82 


20.64 


21.18 


969 


NGC 2997 


Sc(s) L3 


YJHKBrj. 


0.38 


23.8 


31.42 


19.2 


35 


107^ 


32'' 


20.14 


20.33 


21.27 


21.75 


5313 


NGC 4030 


Sbc(r) I 


YJHK 


0.44 


23.6 


31.99 


25.0 


53 


27« 


44« 


19.70 


20.19 


20.76 


21.26 


1196 


NGC 4321 


Sc(s) I 


YJHK 


0.66 


23.7 


32.07 


26.0 


83 


153* 


27* 


19.32 


19.81 


20.74 


21.66 


1184 


NGC 5247 


Sc(s) I-II 


YJHK 


0.41 


23.7 


31.77 


22.6 


45 


97'' 


29* 


19.79 


20.65 


21.26 


21.75 


2259 


NGC 7424 


Sc(s) II.3 


JHK 


0.41 


24.3 


29.88 


9.5 


19 


44' 


27' 


20.80 


21.23 


22.21 




6137 



Notes. References to tiie projection parameters adopted: Zurita et al. (2002), van Zee & Bryant (1999), Lindblad et al. (1997), Pence 
et al. (1990), Hess et al. (2009), <^ Ganda et al. (2006), <«> Knapen et al. (1993), <*> Kuno et al. (2007), and ® Becker et al. (1988). 



shows the strongest correlation with the spiral arms, H^, and ra- 
dio continuum emission, and that this correlation decreases with 
age. 

Most of these studies are based on H„, UV, and visual broad- 
band observations, which are significantly biased by extinction. 
In spiral arms with prominent dust lanes, it is difficult to con- 
duct a complete census of star forming regions and cluster com- 
plexes in the visible especially for the youngest ones, which 
may still be embedded in dust. In these cases, an almost unbi- 
ased, magnitude- limited sample can be obtained more readily in 
the near-infrared (NIR), where the attenuation by dust is much 
smaller (e.g. Ak~ 0.1 Ay). 

In general terms, the K-band light from the disks of spiral 
galaxies is dominated by an old stellar population (Rix & Rieke 
1993). Nevertheless, in a study of the NIR K-band images of 
spiral galaxies Grosb0l & Patsis (1998) noted that several grand- 
design spirals have bright knots along their arms. The alignment 
and concentration of these knots in the arms suggested that they 
were associated with young objects. This was confirmed by the 
NIR spectra of several knots in NGC 2997 (Grosb0l et al. 2006) 
displaying strong Biy emission, which demonstrates that the 
knots are H II regions deeply enshrouded in dust. Both Br^ emis- 
sion and NIR colors can be used as age indicators for sources 
younger than -lOMyr (Leitherer et al. 1999, hereafter SB99) 
and allow the study of the early phases of star formation, which 
are often hidden by dust. 

A sample of 46 spiral galaxies was studied by Grosb0l & 
Dottori (2008), who found that around 70% of these grand- 
design spirals have a significant concentration of bright K-band 
knots in their arms. The NIR color-color diagrams of four of the 
galaxies with JHK photometry suggest that a significant fraction 
of the extended sources are complexes of young stellar clusters 
with ages <10Myr that are reddened by several magnitudes of 
visual extinction. The brightest knots reach Mk - -15™5 corre- 
sponding to stellar clusters or complexes with total masses of up 
to several lO^M©, accounting for the typical SFR of 1 Mq yr"^ 
measured for the brightest grand-design spiral galaxies. 

The main emphasis of our paper is to characterize the general 
distribution of clusters in grand-design spiral galaxies. Relying 
on NIR photometry, this study intends to provide a more com- 
plete sampling than in visual bands including highly reddened, 
young clusters in the arm regions. Our access to younger, partly 



embedded clusters allows us to estimate of extinction at early 
phases of star formation and thereby establish a clearer defini- 
tion of the CLE It is still unclear whether spiral structure can 
influence star formation on small scales (see e.g. Bigiel et al. 
2008; Seigar & James 2002). The combination of NIR colors 
and spatial information offers a detailed view of the distribution 
of clusters as a function of their location relative to spiral arms. 
This is used to test whether strong spiral perturbations can affect 
local star formation. 

A long-term objective, to be pursued in a forthcoming paper, 
is the study of the possible implications of the spatial distribu- 
tion of newly formed clusters on the nature of spiral arms. The 
ability to estimate the individual ages of the youngest clusters 
(Grosb0l & Dottori 2009) could reveal spatial age gradients that 
can be compared with predictions for different spiral arms sce- 
narios. It should be possible to distinguish between density wave 
(Bertin et al. 1989; Gittins & Clarke 2004; Lin & Shu 1964; 
Roberts 1969) and material arms (Wada et al. 2011) owing to 
their different relative angular velocities. In addition, the contrast 
between star formation in arm and inter-arm regions (Seigar & 
James 1998, 2002) could be used as an indicator of the ongoing 
mechanism. However, owing to the intrinsic velocity dispersion 
of clusters and the limited age resolution, the current data cannot 
be used to estimate the longevity of spirals, that is whether they 
are long-lived quasi- stationary spiral structures resulting from 
density waves (Fujii et al. 2011) or transient spiral patterns that 
occur because of recurrent gravitational instabilities (Foyle et al. 
2011;Sellwood 2010,2011). 

In the current paper, we first present the NIR data and de- 
scribe the identification of the stellar cluster complexes and ex- 
traction of their properties. Color-color and color-magnitude di- 
agrams of the cluster complexes are presented in Sect. 4 with a 
discussion of the imphcations for their physical properties. Star 
formation rates and CLFs are derived in Sects. 5 and 6. The gen- 
eral distribution of cluster complexes in the arm and inter-arm 
regions is given in Sect. 7. The appendices includes descriptions, 
maps, and diagrams for the individual galaxies. 

2. Data and reduction 

Eight grand-design galaxies with relatively high SFRs and con- 
centrations of stellar clusters in their arms were selected from 
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Fig. 1. Direct Kg -maps of the eight most prominent grand-design spiral galaxies in the sample. The images are presented in a negative, arbitrary 
intensity scale to more clearly display their spiral structure. All images are orientated with north to the top and east to the left with a scale indicated 
by the 30" bar in the lower right comer. 



the list of Grosb0l & Dottori (2008). The candidates were se- 
lected to have inclination angles lA < 45° and a geometry that 
made it possible to map their spiral structure unambiguously. To 
get a reasonable range of morphological types with strong spiral 
perturbations, including bars, it was necessary to include galax- 
ies with distances up to around 25 Mpc. Two grand-design spi- 
rals with weak spiral arms (i.e. NGC 1232 and NGC 7424) were 
added to allow a comparison of properties as a function of pertur- 
bation strength. The galaxy NGC 7424 has an estimated distance 
of 10 Mpc, which permits a study of its faint clusters and there- 
fore a more reUable estimate of the shape of the cluster distribu- 
tion at low masses. The galaxies are listed in Table 1 where their 
morphological types according to Sandage & Tammann (1981, 
hereafter RSA) are also given. Although only two galaxies are 
classified as strongly barred, all the other galaxies display weak 
bars or oval distortions in their central regions on their K-band 
maps. The eight main grand-design galaxies in the sample are 
shown in Fig. 1 . All the galaxies are displayed in Appendix A, 
where a short description of their morphology is also given. 

The galaxies were observed with the HAWK-I/VLT instru- 
ment, which employs a mosaic of 4 Hawaii 2RG 2048x2048 
chips with a total field of 7' and a pixel size of 0. 106" . All obser- 
vations were performed in service mode except for NGC 7424, 
which was acquired during the HAWK-I commissioning. We de- 
signed a single observing block for each galaxy, including the 4 
broad-band filters, in service mode to ensure that the observing 
conditions were as similar as possible. Target and sky exposures 
were interleaved with a maximum duration of 3™ between empty 
field exposures to allow accurate sky subtraction. Typical total 
exposure times on source were for in the range 9-12™, while 
they were 7™, 4", and 4"" for the filters H, J, and Y, respectively. 
In addition, a Br^-frame with a total exposure of 24™ on source 
was secured for NGC 2997. The individual frames were offset to 
ensure that all parts of the galaxies were exposed including the 
15" gaps between the chips. A minimum of four frames were 



taken for the J- and Y-bands which required shorter exposure 
times. In these cases, parts of the field were covered by only one 
or two exposures, reducing the abiUty to reject bad pixels and 

cosmic ray events there. 

A dedicated, ESO-MIDAS based pipeline was used for the 
data reduction because a generic template, which is not sup- 
ported by the standard ESO pipeline, was used for the ob- 
servation. The observatory-provided dark and flat-field frames 
were used. Since the flats were taken during the morning and 
evening, they may contain small gradients. A low-order, poly- 
nomial correction was applied to the observatory flats estimated 
from stacked sky exposures obtained during the nights. The sky 
frames were inspected visually and rejected if odd features were 
seen. Furthermore, all discrete sources detected on them were 
removed. The average of adjacent sky exposures was then sub- 
tracted from the target frames. The resulting frames had to be 
re-sampled before stacking because the different detector chips 
were not perfectly aligned. A nearest pixel algorithm was ap- 
plied since the pixel size was significantly smaller than the av- 
erage seeing of the frames and an interpolation would smear out 
bad pixels. The transformation constants including gains for the 
individual chips were taken from the instrument manual and ap- 
plied. 

The galaxy frames were aligned to within one pixel using 
stars in the field and centered on the K-band intensity peak of 
their bulges. The quaUty of the individual frames was checked 
by comparing the magnitude and seeing of a bright field star, 
which in a few cases led to the rejection of frames. Bad pix- 
els and outliers could be identified when at least three frames 
were available. These pixels were omitted when their inclusion 
would have significantly increased the variance around the mean 
(Grosb0l et al. 2004). Only for NGC 1566, observed in raflier 
variable conditions, was it necessary to remove a frame from the 
J and Y stacks. This left an unexposed square of 15x15" close 
to the center of NGC 1566 in the J- and Y-maps. 
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Fig. 2. Sources detected in a star forming complex in the southern arm 
of NGC 2997 (a20oo= 9*'45"'33, (Jjooo = -3r21'28") using different de- 
blending contrast thresholds (a) 6c = 0.001 and (b) 6^ = 0.01. Sources 
with identical positions for both 6^ values are indicated by '+', while 
differences are marked by 'A'. 



The HAWK-I field size allowed the use of foreground stars 
from 2MASS (Skrutskie et al. 2006, hereafter 2MASS) for the 
photometric calibration. The uncertainty in the photometric zero 
points was smaller than 07^08 for all maps based on 7-60 2MASS 
stars depending on Galactic latitude of the fields. For the Y-band, 
2MASS magnitudes were estimated from the transformation by 
Hodgkin et al. (2009) for the WFCAM instrument with simi- 
lar filters. Using these 2MASS stars, astrometric transformations 
were also derived with typical errors of <0.2". The full-width 
at half-maximum (fwhm) of stellar images on the final, stacked 
Kj-maps is listed in Table 1 where the surface brightness yUs for 
a signal-to-noise ratio of five per pixel is also given. The other 
bands had similar seeing and surface brightness limits. All non- 
stellar sources were corrected for foreground absorption due to 
our Galaxy using the maps of Schlegel et al. (1998) provided by 
the NASA/IPAC Extragalactic Database (NED). 

The distance moduli for the galaxies were derived from their 
systemic velocities relative to the 3K CMB using a Hubble con- 
stant of 73 km s"' Mpc"' as given by NED. They differ signif- 
icantly in several cases from those obtained using the Galactic 
standard-of-rest (de Vaucouleurs et al. 1991, hereafter RC3), 
such as for NGC 2997 for which the distance adopted places 
the galaxy almost 60% further away. The 3K CMB corrections 
were preferred since they gave more consistent magnitudes for 
the brightest sources in the galaxies. The distances adopted are 
listed in Table 1 together with the linear resolution correspond- 
ing to the fwhm. For most of the galaxies, sources with a size 
of 50 pc could be resolved. This does not allow us to identify 
individual stellar clusters but rather complexes of clusters and 
star forming regions. The scale is compatible with that of typical 
giant molecular cloud (GMC) sizes suggesting that the sources 
detected are associated with individual GMCs. 



T 1 1 r 
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Fig. 3. Instrumental aperture magnitudes for sources detected in the Kj- 
band image of NGC 2997 as a function of their isophotal magnitudes. 
Only stars with photometric errors <0?'2 are plotted. The color indicates 
the value of cs from 1 for stars (red) to for diffuse objects (blue). 



3. Detection of sources 

We performed our search for sources with SExtractor v.2.8.6 
(Berlin & Arnouts 1996) on the K-band frames since they are 
less affected by dust extinction and are therefore expected to 
yield a less biased list of candidates. The steep flux gradients in 
the central parts of the galaxies significantly limited the detec- 
tion of faint sources in these regions. To remove these gradients, 
a smoothed version of the K images was subtracted before apply- 
ing the search. The smoothed image can be generated by using 
either a parametric model fitted to the surface brightness distri- 
bution of the galaxy or digital filtering. The latter option was pre- 
ferred to avoid any bias being caused by the model adopted. The 
filter consisted first of a median filter with a kernel of 5" squared 
followed by a Gaussian smoothing to reduce the noise level. This 
retained all sources with significant gradients on scales smaller 
than the filter size, while removing the smooth variation due to 
the spiral arms, exponential disk, and bulge, except in the central 
5-10" of the galaxies. It allowed the use of a large background 
mesh for SExtractor. The search was performed with a convo- 
lution mask of a fwhm of two pixels corresponding to half of 
the typical seeing. A detection threshold of 2.5 times the back- 
ground noise was used and only sources with at least four pixels 
were accepted. 

Several large and complex sources exist along the spiral 
arms. They were occasionally separated into individual objects 
depending on the SExtractor de-blending contrast threshold pa- 
rameter 6c- A typical complex region in the southern arm of 
NGC 2997 is displayed in Fig. 2, where sources detected with 
two different values of 6c are shown. Without the resolution to re- 
solve individual stellar clusters, the choice of the 'best' value for 
6c is subjective. Comparing several different complex regions, 
6c = 0.001 was found to yield a reasonable separation and was 
used for the final search. 

The classMar estimator, cs, provided by SExtractor was 
used to indicate the shape of sources where cs - \ corresponds 
to star-like images and cs = to diffuse ones. Sources with 
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Fig. 4. Locations of non-stellar sources with errors <(J?3 identified in the galaxies shown in Fig. 1. Compact sources (0.3 < cs < 0.95) are indicated 
by dots (red) while diffuse ones (cs < 0.3) are plotted as triangles (blue). 



cs > 0.95 were assumed to be foreground stars. The distribu- 
tions of non-stellar sources in the eight main targets are shown 
in Fig 4. It is clear that the diffuse sources (cs < 0.3) outline the 
spiral structure more clearly than the compact ones, which are 
more randomly distributed within the galaxies. Only for two of 
the galaxies is this not the case, namely for NGC 1300 where 
the main star formation occurs at the end of its bar, and for 
NGC 7424, which has a weak spiral perturbation. All the galax- 
ies are displayed in Fig. A. 3. 

The centroids of the candidates found by this procedure were 
then used to calculate the magnitudes in each available band. 
Both aperture and isophotal magnitudes were estimated since 
they have somewhat different properties and do not assume a 
specific shape of the source. The former was based on an aper- 
ture with a diameter of 1" and a background area consisting of 
a 2" ring with a radius of 4", to which a /ccr-clipping filter was 
applied. Isophotal magnitudes were computed by mapping the 
Kj-band pixels across a region of 4" around the centroid into ei- 
ther 'source' or 'background' pixels. This was done by moving 
out from the local maximum near the centroid in a spiral pattern, 
while monitoring the pixel values. The background level was es- 
timated as the minimum value within the region for which no 
significant gradient could be determined. Using this level and 
the associated standard deviation, the pixels were classified as 
'source' if they were within the fwhm of or 3cr above the back- 
ground, while all pixels consistent with the level were flagged 
as 'background' . The isophotal magnitude of a source was com- 
puted using this map for all the bands, which ensured that the 
same set of pixels was used to calculate color indices. We in- 
sisted that a source had more 'background' than 'source' pixels. 

The two types of magnitudes are compared in Fig. 3, where 
they are shown for objects detected in NGC 2997. Bright, stellar- 
like sources show a good correlation with a slope of unity. By 
definition, isophotal magnitudes do not record the same fraction 
of flux at all magnitude levels. This effect may approach CTS for 
the faintest magnitude used in this sample. Aperture magnitudes 



may also be brighter when other sources, in crowded fields, are 
included in the aperture. The opposite trend is seen for brighter 
sources, where many non-stellar objects are more luminous ac- 
cording to their isophotal magnitude than that computed for a 
fixed aperture. The effect appears for instrumental magnitudes 
<-8" for which typical non-stellar objects have more 'source' 
pixels than the aperture used (i.e. approx. 280 pixels). This sug- 
gests that the aperture estimate misses flux associated with the 
source because of its fixed size. In general, aperture magnitudes 
give a more reliable and consistent measure of the fluxes asso- 
ciated with the individual sources and were used in the current 
study. The photometric data for the sources detected in the galax- 
ies are available in digital form from CDS'. These tables contain 
the aperture photometry for each filter with associated estimates 
of the error and background. In addition, we provide the coor- 
dinates and cs index. The estimation of the total flux for the 
largest, most complex sources can be obtained more accurately 
from isophotal magnitudes but remains uncertain owing to the 
arbitrary division into substructures (i.e. choice of de-blending 
threshold 6c)- 

For a typical galaxy in the sample with a disk of 10 arcmin^, 
the contribution of background galaxies is 97+12 sources with 
K<20'" or 44+8 with K<19'" using the estimates derived by the 
GalaxyCount program (Ellis & Bland-Hawthorn 2006). To limit 
the number of these galaxies in the sample, only sources within 
the galactic disks were assumed to be defined when their Kj- 
band background surface brightness was 3cr above sky. 

The limiting magnitudes were estimated from the histogram 
of apparent magnitudes as the number weighted average of the 
highest bin and its two neighbors. Comparing with simulated 
fields of sources with a power-law luminosity function, the esti- 
mate corresponds to roughly 90% completeness. The total num- 
ber of non-stellar sources for which aperture photometry could 



' Centre de Donnees astronomiques de Strasbourg: http://cds.u- 
strasbg.fr 
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Fig. 5. Color-color diagram for non-stellar sources in NGC 2997 with 
errors <0';'O5. Evolutionary tracks for a stellar population are shown 
where the full drawn line is for a Padova SSP model while the dashed 
line indicates a SB99 CSP model. The color coding of the points indi- 
cates their absolute magnitude Mk- A standard 'screen' reddening vec- 
tor and one for a 'dusty' cluster medium are shown for Ay = 5™. 



be estimated is given in Table 1, where limiting magnitudes for 
all filters are also listed. 



4. Distribution of magnitudes and colors 

A main tool for the study of the physical properties of the un- 
resolved cluster complexes is their integrated colors and mag- 
nitudes. A typical (H-K)-(J-H) color-color diagram (CCD) is 
shown in Fig. 5 for non-stellar sources in NGC 2997 with pho- 
tometric errors <(F05. For reference, two evolutionary tracks 
are shown for a stellar population with metallicity Z=0.02, 
a Salpeter initial mass function (IMF) (Salpeter 1955), and 
an upper mass limit M„ = 100 M0. The full drawn curve is 
a single-burst stellar population (SSP) based on the Padova 
isochrones (Marigo et al. 2008) including the post thermal-pulse 
asymptotic-giant-branch (TP-AGB) phase. The other track rep- 
resents a continuous star formation population (CSP) from SB99 
indicated by a dashed line. The shift between the two tracks is 
mainly due to the dust and nebular emission that is included in 
the SB99 models. Even though NIR extinction is much smaller 
than in the visual bands, it is still significant. As a first order 
compensation for attenuation by dust, one can construct a red- 
dening 'corrected' color index Q assuming a linear correction 
term. We consider two extreme cases: a screen model of the dust, 
which would yield = (H-K) - 0.563x(J-H) (Indebetouw et al. 
2005), and a mix of stars and dust giving Qd = (H-K) - 0.844x(J- 
H) (Israel et al. 1998; Witt et al. 1992). Both reddening vectors 
are indicated in Fig. 5 for a visual absorption Ay - 5™. 

Two main concentrations of sources can be identified in the 
CCD, namely a) one of the complexes around ((H-K),(J-H)) - 
(0^4, U^S), which corresponds to a stellar population older than 
20 Myr and relatively low extinction, and b) a weaker one close 
to (O'J'8, 1™1), which is compatible with younger and more ob- 
scured complexes. The remaining clusters are more homoge- 
neously scattered to higher (H-K). Their colors can be mainly 



modeled by SB99 adding specific amounts of extinction and 
they are likely younger than 7 Myr. However, to model those 
clusters with the reddest (H-K), one would need to add a sig- 
nificant amount of emission from hot dust. We note that none 
of these young complexes has a visual extinction of less than 
2™. Somewhat depending on the reddening assumed, there is a 
tendency for younger sources to have a higher extinction. 

The color-color diagrams for all the galaxies are given in 
Fig. 6 for comparison. Taking NGC 2997 as a reference, all the 
CCDs share the same basic features, although with different em- 
phasis. The position of the clump at the older end of the cluster 
tracks varies only slightly among the galaxies with the exception 
of NGC 1566. In this galaxy, we recognize the same morphol- 
ogy in the color-color diagram, but all complexes are shifted al- 
most 0^3 towards higher Q values. This may partly result from 
the variable weather conditions under which it was observed as 
mentioned in Sect. 2. The clump of older complexes contains 
a larger fraction of faint sources, which makes its population 
size sensitive to the limiting magnitude. This effect can be ob- 
served by comparing the CCDs of NGC 2997 and NGC 4321, 
where the former has a 1?5 deeper limiting magnitude in Mk. 
There are also intrinsic differences between the galaxies that can 
be seen for galaxies with similar limiting absolute magnitudes, 
e.g. NGC 2997 has more old complexes than NGC 1232, and 
NGC 4321 has more than NGC 4030. 

The concentration of objects with higher (J-K) values, near 
(0^8, 1™2), is present for all galaxies, although the gap to the 
older clump is most clearly seen in NGC 2997, NGC 5247, and 
NGC 7424. Comparing the relative position of the two concen- 
trations with the reddening vectors, it is evident that the objects 
with higher (H-K) cannot be explained by older sources being 
attenuated by dust since they would get higher (J-H) values than 
observed. This shift can best be explained by the sources be- 
ing younger in addition to having high extinction. The young, 
dusty objects in NGC 7424 are shifted even more strongly to 
higher (H-K) values. This could, in part, be due to the higher 
linear resolution, which reduces the probability of several clus- 
ters being included in one target aperture. Without spectroscopic 
data, it is impossible to accurately estimate the physical proper- 
ties of the complexes. Variations in nebular emission and cluster 
IMF can shift the location of sources in the CCD by at least 
0?1. Systematic differences in the amount of bright, red, post 
TP-AGB stars (i.e. 2™ < (J-K)) between the clusters could also 
play a role. Full-scale CCDs for the galaxies are also shown in 
Fig. B.l. 

The evolutionary tracks from SB99 indicate that the 
youngest clusters display strong Br^ and nebular emission. The 
emission decreases with time and vanishes around 7 Myr in a 
way that depends slightly on the model parameters. This sug- 
gests that the Q index, due to the emission in the Kj-band, can 
be used as an age indicator for clusters younger than 7 Myr. This 
was tested using Br^ exposures of NGC 2997, which has a sys- 
temic velocity that permits the HAWK-I Br^-filter to be used. 
The index (Br^-K), with a zero point defined by the foreground 
stars, is shown in Fig. 7 as a function of the Qd index. The rela- 
tion (Br^-K) - -0.48 x Qd - (Fll is determined, with a scatter of 
0?08 for the 55 clusters with significant Bry-emission (i.e. (Br^- 
K) < -0."!) and errors in Bry<0?'03. This confirms that Qd can 
be used to estimate the Br^-emission and the age of clusters with 
O?! < Qd for which the contamination of older clusters is small. 
For Qd<07'l, the index shows little variation and cannot be used 
as an age indicator The equivalent width of Br^ in emission was 
previously calibrated as an age indicator using SB99 models and 
K-band spectra (Grosb0l et al. 2006). 
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Fig. 7. Biy index (Br^-K) as a function of the Qj index for non-stellar 
sources in NGC 2997. Sources with photometric errors in Bry<0'?03 
are displayed with colors indicating Mk- The dashed line represents the 
regression line. 

Color-magnitude diagrams (CMD) for NGC 2997 are given 
in Fig. 8 using the Qd and (J-K) color indices, which are roughly 
orthogonal. For clusters with Qd>(Fl, Qd indicates age, whereas 
(J-K) relates to extinction. An evolutionary track for Padova and 
SB99 SSP models with 10^ Mq mass is shown forreference. The 
exact location of the track depends on both the upper mass limit 
M„ and the upper slope of the IMF, while changes in the metal- 
licity around the solar value are of little importance. The upper 
limit to Qd is closer to (FS, which corresponds well to the SB99 
models with M„ > 60 M0. The relatively few complexes with 
higher Qd values are likely affected by emission from hot dust. 
There is a tail of faint sources reaching negative values of Qd 
around -O'.^S. Many of these are probably globular clusters that 
are similar to Galactic ones, which have an average Qd = -CP4. 
The Qd distribution is narrower for brighter complexes partly 
due to their higher average extinction implied by the CCD. The 
brightest complexes have Mk reaching around -15™5, which cor- 
responds to masses close to 10^ M0 assuming a standard SB99 
SSP model. 




-1 10 12 

Q, (mag) (J-K) (mag) 



Fig. 8. Absolute K-band magnitudes for non-stellar sources in 
NGC 2997 with errors <0P3 as functions of Qj and (J-K) indices. 
Padova and SB99 SSP models for a cluster with IO^Mq are given 
by blue full-drawn and red dashed lines, respectively. The color cod- 
ing indicates the other index, where red corresponds to r?5<(J-K) and 
0';'l<Qd, respectively. The 'dusty' reddening line is shown for Av=5"\ 



The right panel of Fig. 8 displays the (J-K) index, which 
shows more clearly than (H-K) the effects of extinction. For faint 
sources, the lower boundary of (J-K) is sharp, in good agree- 
ment with the tracks for old clusters with Av<l™. The youngest 
sources dominate the high (J-K) values and display a signifi- 
cant scatter due to variable extinction. They form a well-defined 
branch with (J-K) x 1™8 that is clearly separated from the bright, 
older sources with (J-K) around 1™2. The model tracks predict 
that youngest clusters should have (J-K) close to (FO and then 
approach values in the range of 0.5-1™ when massive blue stars 
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Fig. 9. Absolute K-band magnitudes for non-stellar sources in 
NGC 7424 with errors <0'?3 as a function of and (J-K) indices. 
Colors and lines are as in Fig 8. 

Table 2. Brightest member and population size of cluster complexes. 
The fifth brightest absolute magnitude Mk is listed using isophotal mag- 
nitudes, M;;,, and aperture estimates, M^^^. The number of non-stellar 
sources N12 with Mk<-12™ is given with '>' indicating lower limits. 
Exponent a^., and number of sources Nj- of a power-law fit to the lumi- 
nosity function for young and old cluster complexes are listed using the 
subscripts y and o, respectively. 



Galaxy 






N12 


ay 


Nv 


a„ 


N„ 


NGC 157 


-17.8 


-14.6 


453 


1.62 


107 


2.38 


369 


NGC 1232 


-16.5 


-14.6 


354 


2.37 


142 


2.54 


768 


NGC 1300 


-15.9 


-14.8 


214 


1.57 


97 


2.00 


193 


NGC 1365 


-18.7 


-15.1 


705 


1.97 


215 


2.15 


942 


NGC 1566 


-17.8 


-15.3 


>702 


1.75 


339 


2.62 


169 


NGC 2997 


-17.6 


-15.2 


1216 


2.29 


189 


2.49 


1396 


NGC 4030 


-19.0 


-15.4 


>791 


2.25 


45 


2.31 


268 


NGC 4321 


-17.4 


-15.6 


>1066 


1.71 


120 


2.44 


463 


NGC 5247 


-17.6 


-15.1 


974 


2.19 


279 


2.72 


599 


NGC 7424 


-14.1 


-12.6 


40 


1.76 


512 


2.38 


520 



evolve to red super-giants. If the extinction varied only slowly 
with age, one would not expect there to be a separation between 
the young and old complexes in two branches. Thus, the gap sug- 
gests that there has been a rapid reduction of extinction, which 
could happen when the first supernovae explode and dust is ex- 
pelled from the cluster environment (Bastian & Goodwin 2006; 
Goodwin & Bastian 2006; Lada & Lada 2003). 

The CMDs for NGC 7424 are presented in Fig. 9 and are 
similar to those of NGC 2997, except that its cluster population 
is almost 3™ fainter. It also has a group of bright, young clusters 
around (J-K) = 1?8 that is separate from the older ones close to 
1™2. The populations of the two branches as a function of Mk are 
more similar than in the case of NGC 2997. For the faintest mag- 
nitudes below -9™, one starts to see old objects with (J-K)>2", 
which may be individual post TP-AGB stars. The CMDs for all 
galaxies are available in Figs. B.2 and B.3. 



5. Luminosity function of clusters 

Two basic, non-parametric estimates of the luminosity function 
for the population of cluster complexes in the galaxies are the 
brightest complex and the numbers of sources above a given 
magnitude limit. The fifth brightest non-stellar source was used 
as an estimator to reduce the influence of background galaxies 
since fewer than five of these galaxies are expected for Ks<16?5 
according to GalaxyCount. Many of the brightest sources have 
sizes exceeding the 1" aperture used for the standard photome- 
try. Although aperture magnitudes still are the more robust in- 
dicator of the luminosity of compact sources, isophotal magni- 
tudes may give a more realistic estimate of the total flux com- 
ing from the largest star-forming complexes. The fifth bright- 
est complex for each galaxy is listed in Table 2 as M^, and 
M.Ka for the isophotal and aperture magnitudes, respectively. 
Nuclear complexes (i.e. r < 20") were excluded as they are of- 
ten significantly brighter than those in the disk. In NGC 1365 
and NGC 4321 in particular, many complexes in the central re- 
gion were identified with aperture magnitudes reaching -15™5. 
Owing to the inclusion of flux from a larger area, M/^, is 2-3" 
brighter than M^a depending on the morphology of the individ- 
ual regions. As an estimate of the flux originating from a cluster 
complex of a single GMC, the aperture magnitude is accurate 
because of the close match between the aperture used and the 
typical, linear sizes of GMCs. 

The size of the cluster population is estimated from the num- 
ber of non-stellar sources brighter than a given magnitude. A 
limit of Mk=-12™ was used as a compromise between including 
some clusters in NGC 7424 and being brighter than the limit- 
ing magnitude for most of the galaxies. Only three galaxies (i.e. 
NGC 1566, NGC 4030, and NGC 4321) did not fully reach -12" 
at a 90% completeness level (see Table 1). The number of these 
complexes is given as N12 in Table 2. 

The analytical form of the cluster luminosity function is nor- 
mally approximated by either a power law NdL - nL^^dL or 
a Schechter function, which adds an exponential cut-off at high 
luminosities. It is important to avoid bins affected by incomplete- 
ness at the faint end of the distribution as well as high luminosity 
outliers (e.g. background galaxies). The fixed aperture used for 
the photometry can lead to underestimates of the brightest lumi- 
nosities estimated for large complexes (see Fig. 3), which can be 
shifted, for sources with a surface area exceeding the aperture, 
to fainter values that depend on their surface brightness. This 
uncertainty in determining magnitudes for the brightest sources 
made it inappropriate to use a Schechter function to fit the lumi- 
nosity function. 

Even though the K-band is less affected by extinction than 
visual bands, it is still important to correct for dust attenua- 
tion especially for the younger complexes, which may have Ak 
reaching almost 1". The distributions of extinction-corrected ab- 
solute magnitudes Mk„ are shown in Fig. 10, where the cumu- 
lative cluster luminosity functions are indicated by full drawn 
lines. The distributions of the young objects have flatter peaks 
that are shifted to brighter magnitudes than those of the old 
ones: This is mainly due to the higher extinction corrections ap- 
plied to the younger population of complexes. The equal size of 
the two populations for the brightest bins is consistent with the 
mass-independent disruption of clusters found by Chandar et al. 
(2010a) for the Large MageUanic Cloud (LMC). 

The histograms provide evidence of neither a significant high 
luminosity cut-off nor a steepening of CLF slopes for bright clus- 
ters, except possibly for young complexes in NGC 4030 and 
NGC 4321. In the magnitude range where both young and older 
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Fig. 10. Numbers of non-stellar sources as a function of their absolute magnitude Mkq corrected for extinction. Young objects (QJ' 1<Q;]) are plotted 
as red columns while older ones (Qj<0;'l) are shown in yellow. Total numbers of young, Ny, and old. No, objects are written below the name of 
the galaxy. The cumulative cluster luminosity function is shown by a full drawn line. 



cluster populations are complete, they have comparable sizes 
with two exceptions, namely NGC 1566 with an excess of young 
clusters and NGC 4030, which has a deficiency. Power laws were 
fitted to the young and old populations separately with a lumi- 
nosity range estimated visually from the Mk„ distributions. The 
exponents and number of sources used are given in Table 2 and 
show values of a in the range of 1.6-2.4 for the younger pop- 
ulation while the older ones are slightly steeper. It is unclear 
whether this trend is significant because of the corrections for 
extinction, which for the young complexes can reach 1™ in the 
K-band. 

6. Star formation rate 

An estimate of the SFR can be derived from the population of 
young clusters assuming that they follow a given evolutionary 
track. Using a SB99 SSP model with M„ = 100 M©, a Salpeter 
IMF, and metallicity Z=0.02, a cluster with a mass of 1O''M0 
has an absolute magnitude Mk~ -15™6 during first few Myr 
Clusters become fainter by almost a magnitude when their mas- 
sive stars start evolving to the AGB after which they brighten 
again, as seen on the evolutionary tracks in Fig. 8 and Fig. 9. A 
CSP model has a smoother variation in color indices and an in- 
creasing luminosity with time depending on its SFR. The SFR of 
the galaxies, associated with complexes observed above the lim- 
iting magnitude, was calculated by summing up all sources with 
07^1 <Qd assuming that they are younger than 7Myr (as given 
by the CSP model) and their masses are proportional to their lu- 
minosity with 10''Mq corresponding to Mko=-15"5. This value 
gives a good estimate of the cuiTent SFR (i.e. within the last 
7 Myr) since it is little affected by extinction, in contrast to esti- 
mates relying on UV or H„ fluxes which may be absorbed in dust 
lanes. The contribution of clusters below the limiting magnitude 
can be estimated by integrating the CLF power-law fit to the 
young population of complexes down to a given limit. Using the 
exponents listed in Table 2 and a lower cluster mass of 100 Mq, 
the total SFR, was computed. These estimates suffer the uncer- 
tainty given by the power-law extrapolation over several orders 
of magnitude. In contrast, they are independent of the different 
limiting magnitudes and indicate more clearly the relative star- 



formation activity in the galaxies. Both current observed SFR 
and extrapolated total SFR, for the galaxies are listed in Table 3, 
where we also quote SFR estimates for eight galaxies from the 
literature. Our values are two-to-three times higher for six of the 
galaxies, possibly owing to the lower extinction in NIR. A large 
amount of star formation in the central parts of NGC 1566 and 
NGC 4030 may be the reason for our lower SFR estimates for 
these galaxies since crowding and steep gradients could reduce 
the completeness of our sample. The two galaxies also exhibit 
the largest differences between bright complexes as measured 
by isophotal and aperture magnitudes (see Table 2). 



The variation in the SFR per unit area, Ssfr, as a function 
of Galactocentric distance is shown in Fig. 11 using the projec- 
tion parameters listed in Table 1. The eiTors are dominated by 
statistical variations caused by the relatively small number of 
clusters in each bin. The mean, radial Kj surface brightness of 
the disk is indicated using a zero point of 667 Jy (Cohen et al. 
2003). We note that Zsfr follows the K, surface brightness of the 
disk in most cases with the exception of the bar regions (see e.g. 
NGC 1 300, NGC 1 365, NGC 1566, NGC 2997, and NGC 432 1 ). 
The relatively lower Ssfr in the bar regions may be caused by ei- 
ther the radial redistribution of material (e.g. gas) or a lower star 
formation efficiency. Comparing the radial range of the strong, 
two-armed spiral patterns in the galaxies with the radial varia- 
tion in SsFR, one does not see any significant change. This agrees 
with the results obtained by Elmegreen & Elmegreen (1986) and 
Kennicutt (1998) in the sense that a strong spiral perturbation 
does not significantly affect the total Ssfr- 



The surface density of young and older non-stellar sources 
in 10" annular bins is shown in the on-line material Fig. C.l. 
The radial variations in the two populations follow each other 
fairly well in the logarithmic plot indicating that their ratio is 
approximately constant. Only NGC 1566 has a significant excess 
of young sources, which is mainly due to its relatively bright 
limiting magnitudes restricting the number of older complexes 
found. 
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Radius (arcsec) 

Fig. 11. Radial star formation rates for the galaxies in 10" bins with indications of the 10-90 percentile range based on statistical fluctuations. The 
radial range of the main, symmetric spiral pattern is shown by vertical, dashed lines. The average K-band surface brightness of the disk is drawn 
with a dashed-dotted line in units of kjy kpc"^. 

Table 3. Star formation rates and arm/inter-arm populations in the galaxies. The current SFR and the total, extrapolated SFR, for the galaxies are 
given in Mq yr"'. Inner and outer radii (r,, r^) for the region with a symmetric, logarithmic spiral pattern, used for the analysis of arm/inter-arm 
properties, are listed. The percentile values for the Kolmogorov-Smirnov test for arm and inter-arm samples of and Mk being drawn from the 
same distribution are shown as and Pm, respectively. The ten percentile values for the Mk distributions in arm and inter-arms regions are listed 
as M^|(j and M^j^. Numbers of young cluster complexes, N° and N',, in arm and inter-arms regions used to estimate the ratio, Ra ,, of SF in arm to 
inter-arm regions are also listed. Finally, SFR estimates from the literature are shown with references. 



Galaxy 


SFR 


log(SFR,) 


r, 


r„ 


P(2 


Pm 


"'-KlO 


M' 


n;: 


n;. 


Rfl:; 


SFR,,, 


NGC 157 


3.8 


2.16 


20 


56 


0.32 


0.00 


-12.9 


-12.4 


43 


22 


3.5 


0.64" , 0.73" 


NGC 1232 


4.2 


3.61 


22 


36 


0.32 


0.00 


-11.9 


-11.7 


13 


4 


4.5 


2.P 


NGC 1300 


3.6 


1.57 


150 


180 


0.20 


0.93 


-12.8 


-13.0 


16 


10 


1.1 


1.4S 1.7'' 


NGC 1365 


5.5 


2.75 


99 


220 


0.00 


0.00 


-12.8 


-12.7 


129 


37 


2.2 


11.3'' 


NGC 1566 


10.9 


2.86 


49 


99 


0.08 


0.01 


-13.8 


-13.5 


114 


36 


5.6 


9.6^ 


NGC 2997 


10.0 


4.08 


55 


125 


0.37 


0.00 


-12.9 


-12.4 


142 


63 


3.1 


3.7« 


NGC 4030 


1.4 


3.42 


49 


70 


0.45 


0.00 


-13.9 


-13.1 


2 


3 




11.3^ 


NGC 4321 


8.1 


2.62 


99 


134 


0.97 


0.24 


-13.8 


-13.5 


22 


22 


1.2 


5.2'', 18.5^ 


NGC 5247 


13.7 


4.23 


40 


121 


0.16 


0.00 


-13.1 


-12.8 


180 


82 


1.9 




NGC 7424 


0.7 


0.71 


49 


77 


0.70 


0.20 


-10.2 


-10.1 


37 


33 


1.7 





Notes. References to other SFR estimates corrected to current distance scale: Sempere & Rozas (1997), *''* Moustakas & Kennicutt (2006), 
<'■' Omar & Dwarakanath (2005), Schulman et al. (1997), Elmegreen et al. (2009), Thilker et al. (2007), and Kodilkar et al. (201 1). 



7. Spatial distribution 

The positions of sources within the galaxies were computed 
from the projection angles given in Table 1 assuming that they 
were located in the plane of the disk. The shape of the spi- 
ral pattern was derived from ID fast Fourier transforms (FFT) 
of the azimuthal, K-band intensity variation in 1" radial bins. 
The radial range of the main, two-armed spiral pattern was es- 
timated by comparing amplitude and phase of the m=2 compo- 
nent, which is shown in 9-ln(r) diagrams in Fig. A. 2. The regions 
occupied by a symmetric, logarithmic spiral arms are listed in 
Table 3. 

The azimuthal distribution of complexes relative to the two- 
armed spiral is shown in Fig. 12 for NGC 2997, considering only 
the radial range with a strong, symmetric pattern. The upper 
panel displays the absolute magnitude Mk, which for sources 
brighter than -13™ shows a strong concentration in the arm re- 
gions. Fainter clusters have a smoother distribution as a function 
of azimuth that also varies in phase with the arms. A similar be- 



havior can be seen in the lower panel, which shows the Qj index. 
Young complexes with CPl<Qd are mainly located in the arms, 
while the older population has a smaller azimuthal variation. The 
concentration in the arms is stronger for the younger and brighter 
complexes. This general trend is confirmed for all the galaxies 
(see Fig. C.2) but is more clearly seen when the spiral pertur- 
bation is strong, as in NGC 1365, NGC 1566, NGC 2997, and 
NGC 5247. 

To investigate general differences in the cluster distributions 
relative to the spiral pattern, arm regions were defined as a 90° 
azimuthal interval centered on the phase of the m=2 FFT com- 
ponent, yielding equal areas of arm and inter-arm regions. This 
definition allows us to estimate a phase difference among the in- 
dividual clusters relative to the spiral pattern. This is impossible 
when using the relative azimuthal intensity variation to distin- 
guish between arm and inter-arm regions. The contrast between 
arm and inter-arm regions varies significantly, where NGC 2997 
has a larger amount of inter-arm complexes than NGC 1365. 
The arm and inter-arm populations were compared to verify 
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whether the two samples could be assumed to be drawn from 
the same distribution function but merely be of different sizes, 
using the Kolmogorov-Smirnov test and the U-test of Wilcoxon, 
Man and Witney. The distribution function of Qd likely does not 
to change significantly from the arm to the inter-arm regions, 
although the total number of sources varies in phase with the 
spiral perturbation. On the other hand, the distributions of Mk 
are different in arm and inter-arm regions at a 99% level of con- 
fidence. The three exceptions (i.e. NGC 1300, NGC 4321, and 
NGC 7424) all lack signs of strong, coherent star formation in 
the radial regions considered. The Kolmogorov-Smimov statis- 
tics are listed in Table 3 with the U-tests yielding similar results. 
These assumption-free tests are sensitive to all kinds of differ- 
ences in the distributions, such as changes in either the shape 
or mean. Comparing the Mk -histograms, the distributions in the 
arms for the galaxies with strong spiral perturbations show a 
shift towards brighter magnitudes, which is on average 0^4 us- 
ing the ten percentile point of the distributions as an indicator 
(see Table 3). 

We also calculated the SFR in arms relative to inter-arm re- 
gions. The ratios Ra j are listed in Table 3, together with the 
number of young clusters used to calculate this ratio. For the 
galaxies with strong spirals, the SFR in arms is two-to-five times 
higher than in the inter-arm regions, which still have some star- 
formation activity. 

The usage of the m=2 FFT component as a phase refer- 
ence for the spiral is acceptable when considering the general 
arm/inter-arm properties of the cluster distributions. The ac- 
tual arms are imperfectly symmetric and may have azimuthal 
shifts of several tens of degrees relative to their average phase. 
Different models of star formation in the spiral arms (Dobbs 
et al. 2010; Gittins & Clarke 2004; Kim & Ostriker 2002, 
2006) have distinct predictions regarding the phase distribution 
of young clusters relative to the potential minimum of the arms. 
To test these predictions, we need to take into account the phase 
shifts of the individual arms, any asymmetries in the arm pro- 
files, and color gradients caused by age and dust variations. This 
discussion will be postponed to a forthcoming paper dealing 
with the detailed shape of the spiral arms in the galaxies. 

8. Discussion 

The current K-band data are 2-3™ deeper than those presented 
by Grosb0l & Dottori (2008) and have significantly higher res- 
olution, which provide a much fainter detection limit for cluster 
complexes. Only 4 of the 46 galaxies discussed in the 2008 pa- 
per had full JHK photometry. Although the current data are of 
significantly higher quality, it is still of interest to compare them 
with the previous results. The fifth brightest Mk aperture magni- 
tudes agree well, with a spread of 0^' 14, and no significant offset 
after being corrected for the different distance scale used. The 
largest difference of 0^3 was seen for NGC 1566. Some of the 
galaxies (e.g. NGC 1365 and NGC 4321) have bright nuclear 
clusters that are l-2'"more luminous than the ones in the disk. 
The number of non-stellar objects Nn brighter than Mk<-13™ 
was also computed, and corrected for the distance scale. In most 
cases, they were compatible to within 20%, except for NGC 157, 
NGC 1566, and NGC 4030, which all have large, compact star- 
forming regions. The higher resolution of the current data al- 
lows us to distinguish smaller complexes and thereby increase 
the number of detected sources significantly. The SFR estimated 
by Grosb0l & Dottori (2008) was based only on the morphol- 
ogy of the sources, owing to the lack of color information. In 
addition, the number of sources detected was around four times 
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Fig. 12. Absolute magnitude Mk and Qd index of non-stellar sources as 
a function of their azimuthal distance from the two-armed spiral pattern, 
where A© = 0° and 180° denote the relative Kj-band intensity maxima 
given by the phase of m=2 FFT component. The top diagram shows Mk 
with colors indicating the Qd index. The number of clusters with Mk > 
-13™ is displayed as a full drawn curve, while the hatched histogram 
shows the brighter clusters. The low section shows Qd with Mk used 
for color coding. The full drawn histogram indicates older clusters with 
Qd<Wl, whereas younger ones are shown by hatched bars. 

smaller, which explains why the estimates of the SFR are be- 
tween two and four times lower than the current ones. 

The distribution of objects in color-color and color- 
magnitude diagrams is the primary means of understanding the 
physical properties of the cluster complexes. The (H-K)-(J-H) 
diagram allows us to distinguish three populations of sources, 
namely: a) one with colors around (07'4, 0^8) compatible with an 
old stellar population (i.e. >50 Myr) at low extinction, b) another 
with higher (J-K) close to (0'"8, 1™1) indicative of a younger 
population with higher extinction, and c) finally objects scattered 
to higher (H-K) values, the largest part of which can be inter- 
preted in terms of SB99 models plus extinction from dust. To 
reproduce the colors of the objects with the highest (H-K) in this 
group, we would need to add a substantial amount of emission 
from hot dust (see e.g. Witt et al. 1992). The relative importance 
of these three groups varies significantly from galaxy to galaxy. 
The total number of objects in the groups with older sources is 
sensitive to the limiting magnitude because clusters get fainter 
with time. Variations in the current SFR away from the average 
will also influence the relative importance of the old and young 
groups. 
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The detailed interpretation of these distributions depends on 
assumptions such as the stellar evolutionary models adopted, 
nebular emission, attenuation by dust, and the shape of the clus- 
ter IMF including its upper mass limit and slope. The lack of ac- 
curate spectral information for the clusters makes it impossible 
to determine the parameters for the IMF in an independent way. 
We adopted the Padova isochrones (Marigo et al. 2008) as they 
include the late evolutionary phases (e.g. post TP-AGB), which 
are important for the NIR colors. A main difference between the 
two sets of cluster tracks shown in Fig. 5 is either the inclusion 
or exclusion of nebular emission when estimating the NIR col- 
ors for SB99 or Padova tracks, respectively. The inclusion makes 
(H-K) significantly redder in the early phases when hot stars are 
present. The mode of star formation also affects the integrated 
colors, where a continuous formation in the clusters prolongs 
the time-span during which nebular and Br^ emission are sig- 
nificant. Without a good estimate of the intrinsic colors of the 
complexes, the appropriate reddening law can only be selected 
based on general considerations. It is reasonable to assume that 
the attenuation by dust of the integrated light is described by 
a law somewhere between a 'screen' and 'dusty' model (Witt 
et al. 1992) since nearby star forming regions contain significant 
amounts of dust. The reddening vector could even be time de- 
pendent if supemovae are able to remove a major fraction of the 
dust from the inner regions of the clusters. This scenario sug- 
gests that the 'dusty' model should be applied to the youngest 
complexes, while the older ones may be described by a 'screen' 
model. 

The correlation between Bry-emission and the reddening- 
corrected Q index for NGC 2997 (see Fig. 7) clearly indi- 
cates that Q can be used to separate very young complexes (i.e. 
%1 Myr) from older ones. The detailed relation between Q and 
cluster age depends on both the model tracks and reddening vec- 
tor assumed. 

The (J-K)-Mk diagrams show a double-branch structure 
(most clearly separated for NGC 1566, NGC 2997, and 
NGC 7424) for bright complexes with (J-K) » 1?2 and rs, 
respectively. The main branch around (J-K) = l'?2 corresponds 
well to the location predicted for old clusters with relatively low 
extinction. The model tracks suggest that young clusters have 
lower (J-K) values reaching almost O^'O for ages younger than 
one Myr. This region is virtually empty, while the majority of 
the young complexes with 0'?^l<Qd are grouped near (J-K) = 1?8 
suggesting almost 10™ of visual extinction. Since young com- 
plexes must migrate to the branch around 1™2, the gap between 
the two suggests that a rapid reduction in the extinction (i.e. $10 
Myr) occurs just after Br^-emission from the complexes ceases. 
The reddening-corrected Q index shows a more continuous dis- 
tribution, supporting the suggestion that the separation of the two 
branched in (J-K) is due to an expulsion of dust rather than a sud- 
den change in the SFR. 

Information on the early phase of the complexes can be ob- 
tained by comparing the location of the young sources with high 
extinction in the CCD, to the model tracks. The shape of this 
group is mainly elongated along the direction of the redden- 
ing vector. This is inconsistent with a SB99 SSP model and a 
smooth SFR in time, which would yield a much broader distribu- 
tion perpendicular to the reddening vector than observed. Most 
CSP models would be feasible because it would take a consid- 
erable time, several Myr, before the clusters would accumulate 
enough mass to be brighter than the detection limit. During the 
earUest cluster-evolutionary-phase the main sequence lifetime of 
the more massive ionizing stars is clearly comparable to their 
contraction time, which would mimic a continuous star forma- 



tion. Later on, when those stars evolved, Padova SSP models 
represent well the CCDs. Also models with low nebular emis- 
sion would be acceptable in which case the estimated visual ex- 
tinctions could reach Ay = 12"". It is not possible to distinguish 
between these scenarios based on broad band photometry only. 

For comparison, we used the Pessev et al. (2006) sample of 
Small Magellanic Cloud (SMC) and LMC clusters with inte- 
grated colors obtained from 2MASS. Those clusters also show a 
distribution similar to those in our sample of galaxies, indicating 
that our sample also includes intermediate and old globular-like 
objects. Even the peculiar colors, such as those observed in sev- 
eral clusters in the galaxies NGC 7424 and NGC 1365 with (J- 
K)<0P1 and (H-K)>0"5 (lower right part of CCD), can also be 
found in the SMC-LMC clusters, such as NGC 458 (200 Myr), 
NGC 2153 (1.3 Gyr), and NGC 2172 (40 Myr). 

The CLFs can be compared to those derived by Larsen 
(2002) from HST UBV photometry. Their mean exponent for 
the power-law function is = 2.26 + 0.16, while our mean ex- 
ponents for old and young clusters are a,, - 2.40 ± 0.21 and 
ay = 1.95 ± 0.30, respectively. As can be seen, the mean ex- 
ponent derived from the HST observations is more similar to 
that for our old cluster mean exponent. This result is consistent 
with the younger objects having higher extinctions, which would 
probably hide many of them in UV and optical bands analyzed 
by Larsen (2002). We do not find any evidence of the high lu- 
minosity cut-off reported by Gieles et al. (2006). Again here, the 
usage of NIR photometry allows us to obtain a more complete 
census of the youngest and brightest clusters which may be ob- 
scured by dust lanes at shorter wavelengths. Our results are con- 
sistent with Chandar et al. (2010a,b) who estimated an exponent 
of 1.80 for the Magellanic Clouds and found no evidence of a 
high-luminosity cut-off. Since NIR colors provide rather limited 
amount of age information for older sources, we were unable 
to determine whether the properties of the CLF depend on the 
cluster age. 

The distributions of the Mk and Qa of sources as a function 
of Ihe azimuthal phase relative to the spiral arms show a char- 
acteristic variation in the regions of strong, symmetric perturba- 
tions. The total number density of complexes oscillates in phase 
with the K-band surface brightness variations. There is no statis- 
tical evidence for a change in the Qj distribution function from 
the arm to the inter-arm regions. There is clear evidence that the 
Mk distribution in the arms differs from that in-between. The 
brightest complexes are preferentially found in arms. The exis- 
tence of non-linear compressions or shocks in the gas could alter 
the mass spectrum and allow more massive clusters to be formed 
in the arm regions. 

9. Conclusion 

Our analysis of NIR photometry and location of cluster com- 
plexes relative to major arms in 10 nearby, grand-design, spiral 
galaxies had led to the following conclusions: 

- All complexes younger than 7 Myr have significant extinc- 
tion in the range Av=3-7™, while older ones typically have 
Av<l'". 

- The (J-K)-Mk diagrams show two well-separated branches, 
which are indicative of a rapid expulsion of dust around an 
age of 7 Myr. 

- The CLFs are well-fitted by power laws with an exponent 
-1.95 ± 0.3 for clusters younger than 7 Myr and 2.40 ± 0.2 
for older clusters. The distributions do not show evidence of 
a cut-off at their high luminosity end. 
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- Cluster complexes formed in strong spiral arms have a CLF 
that is systematically brighter in Mk by at least (F4 than 
those in the inter- arm regions. 

- The azimuthal average of Ssfr is insignificantly affected by 
strong spiral perturbations, which help to concentrate star 
formation in the arm regions. Outside the grand-design spi- 
ral pattern, star formation is more uniformly distributed in 
azimuth. 

- The SFR in the arm regions of galaxies with strong spiral 
perturbations is 2-5 times higher than in the inter-arm re- 
gions. 

The NIR observations of cluster complexes in grand-design 
galaxies clearly reveal a population of very young clusters highly 
attenuated by dust. These clusters are difficult to access in visual 
bands but are essential for understanding the early evolution of 
clusters. The indication of a shift in their CLF, towards brighter 
clusters, in the arms points to a change in the physical condi- 
tions. This suggests a scenario where star formation is enhanced 
by spiral perturbations (e.g. by increased shear or gas densities, 
or shocks) but in inter-arm regions is due to stochastic processes. 
A full analysis of the formation of clusters in the arms of grand- 
design spirals should include dynamical models for the flow of 
material. 
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Appendix A: Morphology of galaxies and 
distribution of cluster complexes 

The K-band images of the ten grand-design, spiral galaxies in 
the sample are shown in Fig A. 1 as direct images, while 0-ln(r) 
representations are given in Fig A. 2. A short description of their 
spiral structure is given below based on the Fourier analysis of 
the azimuthal intensity variation in their de-projected disks. 

NGC 157: The galaxy has a weak oval distortion at its center, 
which reaches to a radius of around 10" and has a relative 
amplitude of a2 ~ 2%. There is a transition zone to 20", 
where a two-armed, logarithmic spiral with a pitch angle of 
33.4° starts. Its relative amplitude peaks around 40% in the 
range 30"< r <50". The main m=2 pattern terminates close 
to 60", outside which a weaker, tighter pattern exists. 

NGC 1232: It has been argued that this galaxy has a morphol- 
ogy similar to that of our Milky Way (Becker & Fenkart 
1970). It has a bar out to a radius of 16" with 02 = 20%. A 
tight, symmetric, two-armed spiral with (ai, = (0.2, 14.7°) 
emerges from the end of the bar. Its symmetric part termi- 
nates at r = 36", after which the northern arm splits as the 
southem arm continues. A more open pattern can be traced 
to 120". 

NGC 1300: There are indications of a weak central bar with 02 
= 10% within 8". The main bar reaches an amphtude a2 
- 55% around 64", where it widens into ansae. A tightly 
wound set of arms emerges from the end of the bar close 
to 95", while a faint spiral pattem can be traced outside 
150" with (a2, = (0.2, 13.7°). A dynamical study by Patsis 
et al. (2010) suggests that the inner arms are supported by 
stars on chaotic orbits that escape through the Lagrangian 
points at the end of the bar. 

NGC 1365: A spiral structure occupies the central parts, while 
the bar is seen in the range 40-92" with an amplitude 02 up to 
73%. The main spiral starts at the ends of the bar and follows 
a logarithmic spiral with i - 33.1° from 99" to the edge of 
the frame around 220". The amphtude of its m=2 component 
exceeds unity in the outer parts owing to its strongly peaked, 
azimuthal shape. 

NGC 1566: The central part hosts a weak oval distortion with a 
position angle of 30° . Two symmetric regions in the radial 
range of 15-35" and N-S direction show strong star forma- 
tion activity. The arms connect to these regions but the loga- 
rithmic spiral first starts at r = 49" with (02, = (0.7, 27.1°). 
Outside 99", the spiral pattem changes shape. 

NGC 2997: Besides a ring of star forming regions at the center, 
a bar with 02 - 10% ends around 45", where the main spiral 
emerges with (02, = (0.4, 21.2°). The regular logarithmic 
spiral pattern terminates close to 125", where the southem 
arm has an abrupt turn. 

NGC 4030: The inner part has a complex structure with multi- 
ple spiral arms. The main symmetric, two-armed spiral pat- 
tern is located in the region 49-70" with (02, /) = (0.3, 26.5°). 

NGC 4321: The inner part of NGC 4321 is occupied by a cen- 
tral bar inside 13", a main bar to around 67"and several 
ring/spiral structures. A set of prominent spiral arms starts 
at the end of the bar but these are asymmetric (with a much 
stronger southem arm) and tightly wound, in a similar way 
to those of NGC 1300. Close to 99", a logarithmic spiral can 
be traced to 134" with the parameters (02, = (0.5, 29.1°). 

NGC 5247: The spiral pattem of this galaxy, as for that of 
NGC 1566, is a prototype clean, open m=2 logarithmic spi- 
ral. There is evidence of both a central bar within 5" and a 



weak spiral pattem reaching 8" that has the opposite wind- 
ing of the main spiral. After a transition zone, the main spiral 
pattem follows a logarithmic spiral with (a^, i) = (0.6, 34.4°) 
in the range 40" to 121". There is a small phase shift around 
100". The azimuthal profile is peaked with significant am- 
plitudes of higher-order even harmonics. 
NGC 7424: This galaxy is the closest in the sample with a dis- 
tance of 9.5 Mpc. With a linear resolution of 20 pc, most 
sources are likely to be associated with individual clusters. 
There are no traces of either significant dust lanes or large 
star-forming regions. Its bar extends to 38", where a tight 
spiral emerges at an angle of almost 90°. The main set of 
symmetric arms starts around 49" with (02, - (0.3, 12.3°) 
and terminates at 77". Another tight spiral can be traced out 
to 90" with (a2, = (0.1, 9.0°). The spirals have a smooth, 
azimuthal shape with little power at higher frequencies than 
m=2. 

The positions of all non-stellar sources with photometric er- 
rors <0?3 are shown in Fig. A.3, where compact complexes 
(0.3 < cs < 0.95) are plotted as red points and diffuse ones 
(cs < 0.3) as blue triangles. The diflrise objects outhne the 
main spiral stmcture more clearly than the compact ones which 
are more uniformly distributed. This is mainly due to a larger 
fraction of compact sources having Qd<0?l and therefore being 
older on average than the diffuse sources. 

Appendix B: Color and magnitude diagrams 

The (H-K)-(J-H) diagrams of non-stellar sources in the galax- 
ies are given in Fig. B.l, where colors represent their Mk. For 
reference, the 'screen' and 'dusty' reddening vectors and cluster 
evolutionary tracks are indicated in Fig. 5. They have the same 
general features as those described in Sect. 4 but there is a signif- 
icant variation in their relative importance. To a large extent, this 
is due to differences in the limiting magnitude since the 'older' 
clump is dominated by faint sources. The average location of 
the high extinction clump varies mainly in (J-H), which may be 
explained by slight changes in the effective extinction law appli- 
cable to the complexes. 

The Qij-Mk diagrams are presented in Fig.B.2 with the stan- 
dard SB99 SSP model for a cluster of 10^ M©. All the galaxies 
except for NGC 7424 have complexes reaching absolute magni- 
tudes Mk close to -15"". NGC 1365 and NGC 4321 have several 
brighter sources, which all are located in their central regions. 
The brightest complexes are found near Qd=O?'0 but would be 
shifted to sUghtly higher luminosities if extinction corrections 
were applied. The decreasing number of sources moving from 
Qd=0^0 to higher values is caused by the youngest objects, with 
higher Q(j, still forming stars and increasing their luminosity. 
Nebular emission also yields higher Qd measures. Although the 
standard cluster tracks only reach Qd=-0'?3, clusters with values 
up to -ff^l are observed in the LMC-SMC (Pessev et al. 2006), 
which can explain the negative tail of Qd. 

Color-magnitude diagrams for (J-K)-Mk are provided in 
Fig. B.3 with the standard SB99 SSP model as reference. The 
younger population (red/dark) is more prominent for high (J- 
K) values which are indicative of higher extinctions. As dis- 
cussed in the main section, young clusters should have lower 
intrinsic (J-K) values than older ones. That they actually have 
higher values suggests that they are highly obscured. Double 
peaks are clearly seen for NGC 157, NGC 1566, NGC 2997, and 
NGC 7424, while only marginal ones are visible for NGC 1232 
and NGC 5247. A clear gap indicates a rapid expulsion of dust 



p. Grosb0l and H. Dottori: Star formation in spiral galaxies. Online Material p 2 



1 1 I 1 1 1 

NGC0157 


_ T _ 

" NGC1232 " 


^ ■ / 

NGCl30£3^ ^ ^ 

1 i — . 1 ^ .1 ■ — 1 1 1- 


-? \ \ ■ ^ \ 

\ 

\ ' . • 

NGC3.365 / 


- / ' A ^ 
' ^ ' ^ 'i 

I ' 3k ^ \ 

% ' X »■ ■■^ / 

\ ^ ^ - ; 

- . y ■ 

NGC1566 ~ 


NGC2997 

— 1 — 1 — 1 ' .1. — +- — y- 


' - - ^ ^■ 

/ /- "v, V 

/ / ■ ., \ \ 

/ f 
\ ^ _ / - 

H — J " ^" 1 h-'-H r-l 


H— 1 — h |l L . ■ 1 h- 

- ■ 

NGC432r - - 


■■■ ■ 

NGC5247 _i_ 


- . ' ■ - 

x v,gr- iy 

■ T, \ ■ J I- 

NGC 7.424 _ 

1 1 1 ' 1 1 1 1 





Fig. A.l. Kg-maps of all the galaxies in the sample in negative representation. Both intensity and angular scale have been varied to help enhance the 
spiral structure in the galaxies. The two dashed ellipses indicate the radial region of the main, symmetric spiral pattern. All images are orientated 
with north to the top and east to the left with a scale indicated by the 30" bar in the lower right comer. 
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Fig. A.2. 0-ln(r) maps of the galaxies showing the variation relative to the radial Kj intensity profile using the projection angles listed in Table 1 . A 
negative gray scale is used from 0.5 (white) to 1.5 (black). A spherical bulge component was fitted and removed before de-projecting the galaxies. 
The azimuthal angle spans 2n starting from the position angle of the galaxy while the radius begins at 5". The two dashed lines indicate the radial 
range of the symmetric, logarithmic spiral pattern identical to the two ellipses in Fig. A.l. The phase of the m=2 Fourier component is given by 
dots. 
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Fig.A.3. Positions of all non-stellar sources with photometric errors <0';'3 identified on the Kj-maps. Compact sources (0.3 < cs < 0.95) are 
indicated by points (red) while diffuse ones (cs < 0.3) are plotted as triangles (blue). The scale of the images corresponds to that of Fig. A. 1. 
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Fig. B.l. (H-K)-(J-H) diagrams for non-stellar sources with errors <0';'05. The blue full-drawn line show an evolutionary track of a Padova SSP 
model, while the red dashed line indicates a corresponding SB99 CSP model (for details see Fig. 5). The two black lines show the directions of 
the reddening vectors for 'screen' and 'dusty' models (see Fig. 5). Colors indicate the Mk of the clusters. 
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Fig. B.2. Absolute Ks-band magnitude Mk of all non-stellar sources with errors <0?3 as a function of their Qd index. The dashed line shows a 
SB99 SSP model with a mass of 10^ M©. The color code indicates whether the (J-K) index is either above (red/dark) or below 1™5 (green/light). 



p. Grosb0l and H. Dottori: Star formation in spiral galaxies, Online Material p 7 



from young clusters, while the lack of one could be caused by 
both a slower expulsions and a wider spread in attenuation by 
dust. The left edge of the distribution is fairly sharp and corre- 
sponds well to the cluster track. This indicates that many older 
clusters have low internal extinction. 

Appendix C: Spatial distribution of complexes 

The radial surface density of non-stellar sources is displayed in 
Fig.C.l, where red (dark) indicates younger and yellow (light) 
older ones. All galaxies show a steep decline in their outer parts, 
while some (e.g. NGC 157, NGC 1365, NGC 2997, NGC 4321, 
NGC 5247, and NGC 7424) have a relative flat part in flieir disks. 
The distributions of younger and older complexes follow each 
other well indicating that current and past SFRs display Uttle 
radial change. 

The locations of complexes relative to the spiral arms are 
given in Fig. C.2, which plots Mk and Qd against the azimuthal 
distance A6 from the phase of the m=2 FFT component. Only 
sources in the radial range of the strong, grand-design spiral 
pattern (listed in the upper right corner) and photometric errors 
<0^3 are shown. The modulation of the source density with az- 
imuthal angle from the spiral arms is stronger for the galaxies 
with strong, grand-design patterns such as NGC 157, NGC 1365, 
NGC 1566, NGC 2997, and NGC 5247. Both NGC 1232 and 
NGC 7424 have relatively weak spiral perturbations, whereas 
the radial regions for NGC 4030 and NGC 4321 are located in 
their outer parts, which have less star-formation activity. As the 
Kolmogorov-Smimov tests show, the shape of the distribution of 
Qd (i.e. age) does not change significantly from arm to inter-arm 
regions. On the other hand, the excess of bright sources in the 
arms seems real. 
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Fig. B.3. Absolute Kj-band magnitude Mk of all non-stellar sources with errors <0'!'3 as a function of their (J-K) index. The dashed line shows a 
SB99 SSP model with a mass of 10^ Mq. The color code shows whether the Qj index is above (red/dark) or below 1 (green/light). 
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Fig. C.l. Surface density of non-stellar sources with photometric errors <Qf"5. Total numbers of young objects, Ny, (with 0?1<Q) and older ones, 
No, (with Q<Cf?'l) are listed for each galaxy. Red (dark) bars represent young objects, while older ones are indicated in yellow (light). The ratio of 
young to old clusters is plotted as a full drawn curve. 
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Fig. C.2. Absolute magnitude Mk and index of non-stellar sources with errors (r(H-K)<0'^3 as a function of their azimuthal distance Ad from 
the phase of the m=2 FFT component. The radial range used is indicated in the upper right corner and corresponds to the regions occupied by the 
main symmetric, two-armed spiral pattern. Colors from blue to red indicate in the Mk plots the value of Qj, while in the diagrams they show 
Mk as in Fig. 12. 



